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ABSTRACT 

The minimum-energy configuration for the magnetic field above the solar photosphere is curl-free 
(hence, by Ampere’s law, also current-free), so can be represented as the gradient of a scalar poten¬ 
tial. Since magnetic fields are divergence free, this scalar potential obeys Laplace’s equation, given 
an appropriate boundary condition (BC). With measurements of the full magnetic vector at the pho¬ 
tosphere, it is possible to employ either Neumann or Dirichlet BCs there. Historically, the Neumann 
BC was used with available line-of-sight magnetic field measurements, which approximate the radial 
field needed for the Neumann BC. Since each BC fully determines the 3D vector magnetic field, either 
choice will, in general, be inconsistent with some aspect of the observed held on the boundary, due 
to the presence of both currents and noise in the observed held. We present a method to combine 
solutions from both Dirichlet and Neumann BCs to determine a hybrid, “least-squares” potential held, 
which minimizes the integrated square of the residual between the potential and actual helds. We 
also explore weighting the residuals in the ht by spatially uniform measurement uncertainties. This 
has advantages in both not overhtting the radial held used for the Neumann BC, and maximizing 
consistency with the observations. We demonstrate our methods with SDO/HMI vector magnetic 
held observations of AR 11158, and hnd that residual discrepancies between the observed and po¬ 
tential helds are significant, and are consistent with nonzero horizontal photospheric currents. We 
also analyze potential helds for two other active regions observed with two different vector magne¬ 
tographs, and hnd that hybrid potential helds have signihcantly less energy than the Neumann helds 
in every case - by more than 10 32 erg in some cases. This has major implications for estimates of 
free magnetic energy in coronal held models, e.g., non-linear force-free held extrapolations. 

Keywords: Active regions, magnetic helds • Electric currents and current sheets • Magnetic helds, 
models • Sunspots, magnetic helds 


1. INTRODUCTION 

Solar hares and coronal mass ejections (CMEs) are thought to be driven by the release of magnetic energy stored 
in electric currents in the solar corona ( e.g., Forbes 2000 ). Consequently, the origins and structure of coronal electric 
currents are a subject of intense research. Currently, however, measurements of the coronal magnetic held are rare 
and subject to substantial uncertainties {e.g., Lin, Kuhn, and Coulter 2004 ; Tomczyk et al. 2008 ): the vector held as 
a function of three spatial coordinates cannot be directly measured. In contrast, vector magnetograms — maps of the 
magnetic held vector over part of the solar photosphere, typically of active regions — have been made for decades. 
These reveal electric currents normal to the photosphere {e.g., Hagyarcl et al. 1984 ; Leka et al. 1996 ), inferred by 
applying Ampere’s law to the horizontal held (the components tangent to the photosphere). Although they have 
historically been rare, vector magnetograms have recently been made more frequently, thanks to NSO’s SOLIS vector 
spectromagnetograph (VSM; Keller, Harvey, and Giampapa 2003 ), the SpcctroPolarimeter (SP; Lites et al. 2013 ) on 
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the Solar Optical Telescope (SOT; Tsuneta et al. 2008) aboard the Hinode satellite (Kosugi et al. 2007), and the He- 
lioseismic and Magnetic Imager (HMI; Scherrer et al. 2012; Schou et al. 2012) instrument aboard the Solar Dynamics 
Observatory. 

For a given distribution of radial magnetic flux at the photosphere, the hypothetical, curl-free coronal magnetic 
field matching this distribution has the lowest magnetic energy among the set of coronal fields consistent with the 
specified radial field ( e.g ., Priest 2014). Since this minimum-energy field is curl-free, by Ampere’s law it is also current- 
free. Hence, the magnetic energy stored in coronal electric currents is the energy above the minimum-energy state, 
and is referred to as free magnetic energy (e.g., Forbes 2000; Welsch 2006). This energy is available to be released 
impulsively in events like flares and CMEs via dissipation of coronal currents. Increases in free magnetic energy are 
therefore thought to indicate a greater likelihood of flares and CMEs. Since the actual coronal field B throughout the 
coronal volume cannot be directly measured, its energy is usually estimated by some form of modeling. Approaches 
include extrapolation based upon a vector magnetogram (e.g., Dc Rosa et al. 2009) or modeling departures from the 
energy of an initial, current-free state inferred from photospheric magnetic evolution (see, e.g., Kazachenko et al. 2015 
or Cheung and DeRosa 2012). 

Photospheric vector magnetic field observations therefore play a key role in estimating the coronal free energy 
present. Quantifying free energy requires estimating both the energy in the actual coronal magnetic field, B , and the 
energy in the corresponding minimum-energy field. Since the latter is curl-free, it can be represented as the gradient 
of a scalar potential, B p = — Vx- This field is therefore often referred to as a potential field, hence the superscript P. 
The divergence-free condition on magnetic fields implies that x obeys Laplace’s equation, meaning that it is completely 
determined by boundary conditions (BCs) on the potential function. One goal of this paper is to investigate BCs on % 
derived from vector magnetograms. Low (1990) argued that trying to exactly match both the normal- and horizontal- 
field BCs was fundamentally flawed, and suggested using the photospheric normal field as a BC for potential field 
extrapolations. We believe, however, that extra information from observed horizontal fields can be usefully included in 
the potential field solution. Fisher et al. (2010) suggested that solutions which closely match both the normal field and 
the curl-free component of the horizontal field could be used to determine potential field structure in the neighborhood 
of the boundary. We present an approach for determining a potential field that matches both the normal field and the 
curl-free part of the horizontal field in a statistical sense, and can be used to extrapolate a potential field throughout 
the coronal volume. 

In addition, while currents normal to the atmospheric layer imaged in vector magnetograms have been studied, 
few observational constraints have been placed on the presence or structure of horizontal currents within and near 
the photosphere. A second goal of this paper is to quantify differences between the observed and potential magnetic 
fields that indicate the presence of these horizontal currents. Theoretical considerations (e.g., Spruit 1981) imply that 
approximately solenoidal “sheath” currents should be present around the peripheries of active region (AR) flux systems 
in the solar interior, since the coherent magnetic fields that form active regions appear isolated from surrounding plasma 
that typically lacks any spatially coherent field. Sheath currents might also be present at the photosphere. (Although 
strong fields are present in “quiet” regions of the photosphere, their spatial structure is highly intermittent; see, e.g., 
(Sanchez Almeida 2009).) 

The remainder of this paper is organized as follows. In Section 2, we discuss Neumann and Dirichlet BCs for 
extrapolating potential fields. In Section 3, we discuss two approaches for finding potential fields that borrow from both 
BCs: the first weights the normal field and curl-free horizontal field equally; the second, recognizing that uncertainties 
generally differ between the measured field components, uses weighting by parametrized noise estimates to explicitly 
incorporate uncertainties into the field extrapolation. Section 4 explores some implications of our potential held models. 
We conclude with a brief summary of our main results, and a discussion of their implications. 


2. CANONICAL POTENTIAL FIELDS 

On the scale of a typical solar active region, one can approximate the spherical solar photosphere as a Cartesian 
plane. Here, we adopt this approach, defining z to be in the normal direction, with z = 0 at the photosphere. In 
the three-dimensional (3D), Cartesian, half-space satisfying z > 0, the potential held B p = — Vx, where the scalar 
potential x satisfies 


x(x,y,z) 
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solves Laplace’s equation if k z = yJk x + k y = kh, with the outer BC set by |£? p | —> 0 at infinity. This functional 

form implies that the spectral function x(k) is a function of only the two independent wave numbers k x and k y : the 
variation with height z is determined by the spatial variations of the potential function in the z = 0 plane. Typically 
either Neumann or Dirichlet BCs are imposed at z = 0 to determine x{k). 

For the Neumann condition, the solution for the potential xn on the z = 0 boundary is then found from 


B z (x,y,0) = -0 z xnU=o = 


+ oo 


/»+oo 
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d k x 


dk v Mkx,k v )k h e ik * x+ik y y 
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Xn (k x ,k y ) = 
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(3) 


where we have defined B z {k x ,k y ) as the Fourier transform of B z (x,y). The spectral function determined by Equation 
(3) can then be used in Equation (1) to determine the potential held for z > 0 that satisfies the Neumann BC. 

We now turn to the Dirichlet BC. In principle, knowledge of just one component of the measured horizontal held 
could be used to specify a Dirichlet BC: the curl-free condition on Bh implies k x B y = k y B x , where B y and B x are 
the Fourier transforms of B y and B x , respectively. This could be used to determine the 2D potential function on 
z = 0, a Dirichlet condition for the 3D Laplace’s equation on z > 0 (see below). Such an approach would, however, 
ignore observations of the other horizontal component, and would therefore generally be inconsistent with these. In 
addition, this approach would also treat the data as exact, without accounting for the presence of uncertainties in the 
measurements. 

Information from both components of the horizontal held are combined in computing the horizontal divergence of 
B at the photosphere, which gives 

V fc • B h = -(d 2 xX n + d 2 X v) ■ (4) 

This is a 2D Poisson equation for xd(x, j/) at z = 0. Specifying Dirichlet or Neumann BCs on xd for the x and y 
boundaries (possibly at inhnity) uniquely determines xd on z = 0, to within a constant (Jackson 1975). Since B p 
only depends on derivatives of the potential function, this forms a Dirichlet condition for the 3D Laplace’s equation 
for xd for z > 0. The solution xd on z = 0 is then found from 


V,> • B h (x,y,0) = -{d 2 XD + <9yXD)U=o 
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where i kh ■ Bh{k x ,k y ) is the Fourier transform of • Bh{x,y, 0). (The divergence-free condition on magnetic helds 
implies that \Bh(x, y, 0)| decays sufficiently fast with x and y that this integral converges for a localized source, such 
as an active region.) This spectral function xd on z = 0 can then be used with Equation (l) to determine the potential 
field for z > 0 that satisfies the Dirichlet BC. 

Vector magnetograms can therefore be used to determine spectral functions for the potential field in at least two 
ways: using the Neumann BC to derive %n from B z , and using the Dirichlet BC to derive xd from V/, • Bh, using 
Equations (3) and (7), respectively. If the magnetic field on the z = 0 boundary is potential, then the uniqueness 
of solutions to Laplace’s equation (Jackson 1975) implies that Vxd (k x ,k y ) = Vxn (k x ,k y ) and therefore B p and 
V/j • B p = —d z B p are consistent. Even if the observed field B ohs were potential, however, the presence of noise or 
systematic errors in the data would, in general, introduce inconsistencies between xn and xd for those observations. If, 
in addition, B ohs is not potential, then generally no potential field will be simultaneously consistent with V/, • B j( bs and 
B° hs . Differences between xn an d xd might be useful to infer properties of either (i) noise / errors in the measurements 
or (ii) currents, or both. 

We note that although our expressions for the Neumann and Dirichlet solutions are only appropriate for a particular 
domain (the Cartesian, z > 0 half space), a different choice of domain would not affect the general conclusion: 
the observed normal field B n and V/, • Bh at the photosphere separately determine Neumann and Dirichlet BCs, 
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respectively, for a solution y to Laplace’s equation in the coronal domain, and the potential magnetic fields derived 
from each of these BCs will, in general, differ due to measurement errors and non-potentiality at the surface. 

We implemented codes to rapidly compute a potential field from either Neumann or Dirichlet boundary conditions on 
Cartesian domains using discrete Fourier transforms (DFTs). We verified the accuracy of our reconstructed potential 
field (and its associated magnetic energy) using synthetic photospheric magnetic held data from a test case with a 
relatively easy analytic solution: the potential held from a submerged, vertically-oriented dipole. 

To demonstrate that differences in the Neumann and Dirichlet potential helds are present in solar data, we solve 
Laplace’s equation for each BC for a SDO/HMI vector magnetogram of AR 11158 from the sequence analyzed by 
Welsch, Fisher, and Sun (2013). We arbitrarily chose an observation from 15 February 2011, with the center of the 
integration time at 02:00 UT (near the end of an X2.2 hare), when the center of the active region was at S20 W13. In 
images of the active region’s helds, we show a 426 x 276 pixel region containing essentially all of the region’s flux; the 
full magnetogram is 612 x 610 pixels, and the area outside of the central region lacks any signiheant spatially coherent 
magnetic held. These data have been interpolated onto a Cartesian plane using a Mercator projection (Welsch et al. 
2009), with a pixel size of 362 km. The ratio of total signed hux to the total unsigned flux over the magnetogram is 
< 1 %. 

Figure 1 compares the observed B z and V;, • Bh, and the versions of these reconstructed from the Dirichlet and 
Neumann solutions, respectively, as well as the difference between observed and reconstructed quantities. Note that 
the saturation levels in this hgure differ between left and right columns, and are smaller in the bottom row. We have 
expressed the inverse length scale in V/,. ■ Bh as pix -1 to facilitate comparisons with the magnitude of B z . Note the 
opposite-polarity halos present around strong-held regions in the middle-left panel, which shows B z derived from the 
Dirichlet BC. This systematic feature of Dirichlet-derived normal potential helds cannot plausibly have arisen from 
random errors in the magnetogram measurements, although they might be due to systematic errors in the inversion 
process (discussed further below). We also verihed that these halos are not an artifact of periodicity inherent in our 
use of DFTs to solve for the potentials, by padding the 621 x 610 pixel magnetogram arrays with a moat of 400 
zeros. (The 426 x 276 arrays shown in this and other hgures were cropped close to the active region to display its 
helds more clearly.) The vertical components of the Dirichlet potential helds derived from the unpadded and padded 
magnetograms closely match: the mean and maximum unsigned difference between the two are 2 and 7 Mx cm -2 
(Mx: maxwell), respectively, over the 426 x 276 held of view shown here. Further, the halos, which are coherent over 
tens of pixels, do not exhibit the high spatial frequencies of Fourier ringing / Gibbs phenomena. (See also Figure Al 
in Appendix A, which shows halos also present in B z in a Dirichlet potential held computed using a non-periodic, 
finite-element approach.) Assuming the l? 2 -halos reflect a real property of V/, • Bh, they indicate that the observed 
vertical held is not consistent with the Dirichlet potential held — i.e., the lack of halos in the observed B z indicates 
that horizontal currents are present (also discussed further below). 

3. HYBRID POTENTIALS: USING BOTH B z AND V h ■ B H 

As a starting point, we assume that the potential held should be the current-free held that is, in some sense, most 
consistent with observations. As we have seen, the Neumann and Dirichlet solutions are each consistent with distinct 
components of the observed photospheric magnetic held. We now consider generalized dehnitions of “most consistent 
with observations.” 


3.1. Least-Squares Difference from the Data 

Instead of computing a potential held that matches either the vertical or horizontal parts of the observed pho¬ 
tospheric magnetic held, another approach is to determine a single potential held that is most consistent with all 
three components of the observed magnetic vector. Adopting the sum of squared residuals between the photospheric 
potential and magnetogram vector helds as our measure of such consistency, we can determine a “combined” potential 
Xc whose gradient minimizes the corresponding functional I, 

/ +oo r+oo 2 

dz / d y ( B° hs + V Xc ) , (8) 

-oo J — oo ' ' 

where all quantities are evaluated at z = 0. As posed, this is a type of least-squares minimization problem, albeit 
non-local. 

Given a variation y 1 about the extremal potential \ c that minimizes I, where y 1 also obeys Equation (1), the 
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Figure 1. Top row: Observed B z (left) was used as a Neumann condition to derive Vh • Bh expected from a potential field 
(right). Middle row: B z (left) for the potential field inferred using the observed X7 h ■ Bh (right) as a Dirichlet BC. Bottom left: 
Difference between observed and reconstructed B z (top left minus middle left). Bottom right: Difference between observed and 
reconstructed Vh ■ Bh (middle right minus top right). The x and y axes are in interpolated, reprojected HMI pixels, which 
are square and 362 km on a side. In the left column, image saturation is (top to bottom): ±1000 Mx cm -2 , ±1000 Mx cm -2 , 
and ±500 Mx cm -2 . In the right column, image saturation is (top to bottom): ±333 Mx cm -2 pix -1 , ±333 Mx cm -2 pix -1 , 
and ±167 Mx cm -2 pix -1 . Note the opposite-polarity halos around strong-field regions in the image of B z derived from the 
Dirichlet condition (middle-left panel). 


first-order variation SI must vanish for arbitrary X '. The first-order variation is given by 
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In the final expression, h points horizontally outward along the closed curve defining the integral. We assume this 
integral vanishes due to either the finite spatial extent of the observed field’s horizontal divergence or constraints on 
the imposed variation x! at the x and y boundaries. 

We Fourier transform to a function of k! h , writing its spectral function as x', giving 


p+oo 


p+oo 


p+oo 
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dx 


d y 


d k’ 


/* + 00 


d k'e ik > +ik 'y y x 


-X'(k' h )k' h {B? s + d zX c)~x'(K)( 
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V h - B° bs + V 2 hX c) 


( 12 ) 


where the differential operator d z acting on x!, evaluated at 2 = 0, has brought down a factor of —k' z = —k' h . Next, 
we Fourier transform B ° hs , Xc , and B^ s into functions of kh, as above using tildes to denote spectral functions, and 
then act on these with the differential operators d z , V^, and V 2 , which bring down factors of —kh, i kh, and — k 2 , 
respectively, to give 


-i /* + oo /*+oo /* + oo /* + oo /*+oo /*+oo 

51= —j / dx dy dk' x dk' y e ik * x+ik v v / dk x / dk y e ik * x+ik y y x 
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We then interchange the spatial and wavenumber integrations, and note that 
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and similarly for the dy integral. Hence, integrating over x and y implies k' h = —kh, so 


1 r+oo r+oo f , . , , 

= 2^2 J dkx J dk y x'{-k h ) k h {B° ha (k h ) - k h Xc{k h )j+ik h ■ B° h s (k h ) - khXc(k h ) ■ (15) 


Since 51 must vanish for an arbitrary variation x' about the extremal potential y c , the quantity in square brackets 
must be zero, i.e., 

0 = k h j. B° hs (k h ) - k h x c (k h )| + i k h ■ B°^ s (k h ) - kl Xc (k h ) , (16) 

which implies 

Xc(fc/,) = [k h B° hs (k h ) + i k h ■ Bt\k h )]/{2kl) . (17) 

Hence, this spectral function can be determined directly from observed quantities (i.e., the spectral functions of B° hs 
and B )( bs ), from which this combined potential field can be calculated. Comparison of Equation (17) with Equations 
(3) and (7) shows that the combined result is an average of the Neumann and Dirichlet spectral functions. This result 
might be expected, from the known least-squares properties of averaging, but it nicely confirms our earlier argument 
that the horizontal divergence is a good choice for specifying the Dirichlet BC, as in Equation (6). 

We have used this combined spectral function with DFTs to compute the combined potential field for the vector 
magnetogram from Figure 1. To confirm that the combined field does agree more closely than either the Neumann 
or Dirichlet fields, we compare the sums of squared differences between each potential field and the observed field, by 
substituting each potential function in for x in Equation (8). We find /n = 1.74 x 10 25 Mx 2 cm -2 , Id = 1.77 x 10 25 
Mx 2 cm' 2 , and I c = 1.36 x 10 25 Mx 2 cm -2 , for the Neumann, Dirichlet, and combined fields, respectively. In Figure 
2, we show the observed B z and V/i • Bh in the top row, reconstructions of these from the combined potential solution 
in the middle row, and the difference between observed and reconstructed quantities in the bottom row. Note that 
the image saturation in the bottom row of Figure 2 is set lower than that in Figure 1, reflecting the fact that the 
combined solution agrees more closely with the observations. In this figure, note that the opposite-polarity halos 
around strong-field regions in the image of B z from the combined field (middle-left panel) are discernible, but weaker 
than those the middle-left panel of Figure 1. We will discuss the physical significance of differences between potential 
and observed fields in more detail in Section below. 

In the top panel of Figure 3, we show histograms of differences in (signed) B z between the observations and each 
of these three potential fields — Neumann (red), Dirichlet (blue), and combined (black solid). Differences are small 
in the Neumann case, since the Neumann field is constructed to match B z ; its histogram has been rescaled to the 
maximum of the combined field’s. Differences are large, however, in the Dirichlet case. As expected, the combined 
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Figure 2. Top row: Observed B z (left) and Vh ■ Bh (right). Middle row: Combined potential field’s B z (left) and Vh ■ Bn 
(right), derived from the spectral function in Equation (17). Bottom left: Difference between observed and reconstructed B z 
(top left minus middle left). Bottom right: Difference between observed and reconstructed Vh ■ Bh (top right minus middle 
right). The x and y axes are in interpolated, reprojected HMI pixels, which are square and 362 km on a side. In the left column, 
image saturation is (top to bottom): ±1000 Mx cm -2 , ±1000 Mx cm -2 , and ±250 Mx cm~ 2 . In the right column, image 
saturation is (top to bottom): ±333 Mx cm -2 pix -1 , ±333 Mx cm -2 pix -1 , and ±83 Mx cm -2 pix -1 . 


field’s differences fall between the other fields’ differences. A Gaussian fit to the combined distribution is plotted with 
a black dashed line, and its fitted width is printed. This can be used to relate the magnitude of field differences to 
noise levels in the observed quantities. In the bottom panel, we show analogous histograms of differences in • Bh 
between the observations and each of the three potential fields — as above, Neumann is in red, Dirichlet in blue, and 
combined in solid black. In contrast to differences in B z , now it is the Dirichlet case in which discrepancies are small 
(since the Dirichlet field is constructed to match V/ t • Bh), and its histogram has been rescaled to the maximum of 
the combined field’s. As expected, differences from the Neumann case are large, and differences from the combined 
field fall between the two other two cases. Again, a Gaussian fit to differences from the combined field is overplotted 
with a black dashed line, and the fitted width is printed. (Our use of DFTs to solve for Dirichlet and Neumann fields 
yields flux-balanced fields to about 1 part in 10 6 or better.) 

Equation (17) determines the spectral function of the potential field that satisfies the particular definition of “most 
consistent with observations” — the minimized, integrated, squared discrepancy — given in Equation (8). Since any 
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Observed — Potential, B z [Mx cm 2 ] 



Observed — Potential, (V h - B h ) [Mx cm 2 pix '] 


Figure 3. Top: Histograms of observed minus potential B z , for Neumann (red), Dirichlet (blue), and combined (black solid) 
potential fields. The Neumann field’s histogram has been rescaled to the maximum of the combined field’s. Bottom: Analogous 
histograms of observed minus potential Vh • Bh for Neumann (red), Dirichlet (blue), and combined (solid black) fields. Here, 
the Dirichlet field’s histogram has been rescaled to the maximum of the combined field’s. As expected, the combined field’s 
difference histograms fall between the other fields’ histograms for each variable. In both cases, a Gaussian fit to the combined 
distribution is plotted with a black dashed line, and its fitted width is printed. The width is a measure of the significance of 
these field differences relative to noise levels in B z and Vh • Bh; see text. 

other potential function (other than a trivial constant function) added to x'c will alter the total potential and/or its 
normal derivative on the boundary (Jackson 1975), this least-squares potential is unique. As we have seen, either 
the Neumann or Dirichlet potential can closely match either vertical or horizontal components of the photospheric 
field, respectively. But in a statistical, least-squares sense, neither matches the total photospheric field, over the full 
magnetogram, as closely as the combined solution. Also, even where they do match the observations, the Neumann 
and Dirichlet approaches have a clear drawback in that they overfit the data: the noise certainly present in the 
fitted component of the observed field is matched precisely (in principle, to machine round-off error), as though it 
were real signal. This is an issue of precision versus accuracy: either the Neumann or Dirichlet solution can match 
the corresponding field component more precisely than the accuracy of its measurement. In contrast, the combined 
potential field matches the observations statistically. 

3.2. Uncertainty-Weighted, Potential Fields 

Uncertainties in measurements of line-of-sight (LOS) fields are typically smaller than uncertainties in transverse 
fields ( e.g ., Hoeksema et al. 2014; Kazachenko et al. 2015). Consequently, for magnetograms of active regions near 
disk center, where the LOS and photospheric normal directions coincide, one could argue that the vertical photospheric 
field is measured more accurately than the horizontal field, and, therefore, that the Neumann BC should be preferred 
over the Dirichlet condition. 
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Two counterarguments can be made. First, even if uncertainties near disk center are much smaller for B z than 
for V;, • B} j, completely ignoring V/, ■ Bh is unwarranted: even if noisier, the horizontal field measurements do 
contain additional information about the structure of the photospheric magnetic field. Given that the observations 
are uncertain, overdetermining the solution by using all available data should yield a result more consistent with 
observations. Second, uncertainties in the Neumann condition increase away from disk center, since the measured 
normal field includes both LOS and transverse components. Hence, noise levels for the Neumann condition can 
become comparable to those for the Dirichlet condition. (We note, however, that deriving the Dirichlet condition 
involves differentiating the data, which will increase uncertainties in the BC: a four-point finite-difference stencil 
for the horizontal divergence effectively doubles the error in the horizontal divergence compared to the error in the 
horizontal field.) 

These considerations raise the question: in the presence of differing noise levels in the vertical and horizontal fields, 
how should these data be used when computing the potential field? Intuitively, we expect the BC should involve both 
vertical and horizontal fields, but these should be weighted by their noise levels. Following the approach outlined 
above, we seek a “noise-weighted potential” Xw whose gradient minimizes a functional / w equal to the integral of the 
weighted sum of the squared residuals between the potential field and the observed field, with differing weights w z and 
Wh for B z and Bh , respectively, 


C+oo 


r»+oo 


iw = 


dx 


d y w z (B° bs + d zXw Y + w h (B° bs + V, lXw ) 


(18) 


One can then seek conditions for which the first order variation 6I W vanishes for any variation X ' about the noise- 
weighted potential Xw . 

One might use weights w z and Wh that vary in space, for instance to account for either geometric dependence of 
uncertainties in the LOS and transverse measurements as functions of disk position or uncertainties propagated through 
the field inversion process in each pixel. Potential magnetic fields are, however, non-local, so errors in the underlying 
field measurements will necessarily be propagated, perhaps globally: a poor measurement of vertical flux density in a 
pixel will affect the potential solution outside of that pixel, even if that particular pixel’s field is not strongly weighted 
in the fit. As a practical matter, use of spatially varying weights complicates the analytic minimization of / w . While 
numerical approaches might enable determination of a potential field that minimizes / w with spatially varying weights, 
for simplicity we assume w z and Wh are constants here. The analysis leading to Equation (16) is then essentially 
unchanged, and leads to the constraint equation on the spectral function Xw (kh) 


0 = w z k h 


B? s (k h ) - k hXw {k h ) 


+ w h 


i kh ' Bh (kh) kh X w[kh) 


(19) 


which implies 


If we set 


obs . 


Xw(^/l) — 


[w z k h B° bs (k h ) + w h ik h ■ B h (k h )\ 


Wh 


k 2 

K h 


1 

W z = — 
<J Z 

1 

Wh = - , 

CTh 


( 20 ) 


( 21 ) 

( 22 ) 


where a z and ah represent the uncertainties in B z and As(V/, • Bh ), respectively, and As is the width of pixels 
(assumed square), then 


Xw(fcft) 


obs . 


[a h k h B° bs {k h ) + a z \k h • B h (k h )\ 


<?h 


K 


(23) 


To explore this weighting with real data, we estimated noise levels from histograms of B z and • Bh from the 
SDO/HMI magnetogram, assuming the cores of each distribution arise from noise. We defined the cores to be values 
within ± 20 Mx cm -2 and ± 60 Mx cm -2 pix -1 , respectively, and fitted each core to a Gaussian. Figure 4 shows the 
distributions and fits. The fitted widths for B z and V/j • Bh are 15.8 Mx cm -2 and 61.6 Mx cm -2 pix -1 , respectively, 
which we interpret as the noise levels in each quantity. (Varying the widths of the fitted core regions did not change 
the fitted widths significantly.) This implies weights of about 4:1 on the Neumann versus Dirichlet condition. 
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-200 -100 0 100 200 
B z [Mx cm -2 , red] & (V h - B h ) [Mx cm -2 pix -1 , blue] 


Figure 4. Histograms of observed B z (red solid) and Vh ■ Bh (blue solid). Gaussian fits to the core of each distribution (within 
± 20 Mx cm -2 and ± 60 Mx cm -2 pix -1 , respectively) are overplotted as dashed lines in corresponding colors, and the widths 
of each fit are given. At this time, AR 11158 was at S20 W13. 


We then used these noise levels in Equation (23) to derive a noise-weighted potential % w and its associated fields. 
Evaluating the sum of squared differences in Equation (8) using % w yields / w = 1.56 x 10 25 Mx 2 cm -2 , so the weighted 
potential agrees more closely with the observed field than either the Neumann or Dirichlet fields, but does not agree as 
closely as the unweighted combined field. (As noted above, In = 1-74 x 10 25 Mx 2 cm -2 , Id = 1.77 x 10 25 Mx 2 cm -2 , 
and I c = 1.36 x 10 25 Mx 2 cm -2 .) This field’s B z does not look significantly different than those from the combined 
field shown in Figure 2, so we have not included a separate figure to show it. Histograms of the residuals, however, 
are different, as we show in Figure 5. It can be seen that the relative weighting has, as expected, improved consistency 
with the observed B z , at the expense of decreased consistency with the observed S7h • Bh- For each observable, the 
fitted widths of the cores of the distributions (±20 Mx cm -2 for B z , ±60 Mx cm -2 pix -1 for • Bh) of the residuals 
are similar to our estimated noise levels. We note that the weights of each observable could be changed to vary the 
widths of the residuals’ distributions. Since typical weighted least-squares fits usually weight by the square of the 
uncertainties in the fitted variable, we also tried this weighting. Given the roughly 1:4 ratio of noise levels between 
the vertical and horizontal fields, respectively, this gives a weighting of 1:16, heavily favoring the Neumann solution. 
With this “squared noise weighting” (snw), the sum of squared differences from Equation (8) yields / snw = 1.68 x 10 25 
Mx 2 cm -2 , worse overall than weighting by the first power of each uncertainty. 

How do the extrapolated fields differ, qualitatively? A selection of field lines is plotted in Figure 6, for both 
the Neumann and noise-weighted hybrid potential fields. In the figure, field line integrations were initialized from 
the same footpoints with the same colors, but only field lines with both footpoints anchored at the photosphere are 
plotted. Many field lines look essentially the same, but some differences can be discerned. We note that the periodicity 
implicit in the DFT method used to compute these potential fields can affect field line connectivities. Minor differences 
aside, the fields’ structures appear similar. Since Neumann potential coronal fields in some cases agree qualitatively 
with observed coronal loops (see, e.g., Schrijver et al. 2005), it is reassuring that the noise-weighted coronal field 
exhibits similar morphology. We remark that observations of coronal magnetic field structure cannot readily be used 
to determine whether one approach to potential field extrapolation is “better” than another, since the coronal field 
is presumed to not be potential. We defer detailed investigation of differences in the fields’ connectivities to a future 
study. 

We have minimized a weighted sum of squared residuals, in a manner analogous to the weighting used in least- 
squares fitting. As with the unweighted minimization, for a fixed choice of weights, the solution to this minimization 
is unique. We separately weighted the horizontal- and vertical-field residuals by our distinct, constant estimates of the 
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Observed — Potential, B z [Mx cm 2 ] 



Observed — Potential, (V h * B h ) [Mx cm z pix '] 


Figure 5. Top: Histograms of observed minus potential B z , for Neumann (red), Dirichlet (blue), and noise-weighted (black 
solid) potential fields. The Neumann field’s histogram has been rescaled to the maximum of the combined field’s. Bottom: 
Analogous histograms of observed minus potential S7h ■ Bh for Neumann (red), Dirichlet (blue), and noise-weighted (solid 
black) fields. Here, the Dirichlet field’s histogram has been rescaled to the maximum of the combined field’s. A Gaussian fit to 
each noise-weighted distribution is plotted with a black dashed line, and its fitted width is printed. The widths of the noise- 
weighted distributions are commensurate with the noise levels in the observed quantities, shown in Figure , used to determine 
the BC. 

uncertainties in each. 

4. PHYSICAL SIGNIFICANCE OF HYBRID BOUNDARY CONDITIONS 

4.1. Implications for Free Energy Estimates 

The free magnetic energy in the coronal magnetic field is the difference between the energy in the actual coronal 
field and the energy of a hypothetical, minimum-energy field. It can be shown that the Neumann potential field is the 
lowest-energy field consistent with a given normal magnetic field at the photosphere {e.g., Priest 2014). 

While B cannot currently be measured throughout the coronal volume above an active region, departures from the 
potential state can be estimated by either modeling B in the corona {e.g., De Rosa et al. 2009; Cheung and DeRosa 
2012), or quantifying the photospheric Poynting flux {e.g., Kazachenko et al. 2015). Quantifying coronal free energy 
is a key goal of such efforts. 

Do different choices of BCs for potential fields result in different energies for the potential field? Simply put: yes. 
The energy U of a potential field can be written as a surface integral of the potential function times the normal 
magnetic field {e.g., Welsch 2006) at the photosphere, 

u= hj iAB ’- x ' 


(24) 
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Figure 6. Top: selected field lines integrated in the Neumann potential field. Bottom: selected field lines integrated in the 
noise-weighted potential field. Integrations were initiated from the same footpoints, and are plotted with the same color. The x 
and y axes are in interpolated, reprojected HMI pixels, which are square and 362 km on a side. All field lines shown have both 
endpoints at the base of box. 


where B z is the normal field that is consistent with the potential function %, so B z in Equation (24) differs from 
B° hs for all potentials but the Neumann solution. Using the observed vector magnetic fields shown in the top row 
of Figure 1, the energies of the Dirichlet, Neumann, combined, and noise-weighted potential fields (derived using the 
spectral functions in equations 3, 7, 17, and 23, respectively) are 4.6 x 10 32 ,9.1 x 10 32 ,6.1 x 10 32 , and 7.7 x 10 32 
erg, respectively. Because the Dirichlet field’s energy is roughly a factor of two lower than the Neumann field’s, the 
combined and noise-weighted potential-field energies are also significantly smaller than the Neumann potential-field 
energy. It seems plausible that the excess energy in the Neumann solution arises because the observed normal field 
itself is, in a key sense, non-potential. This is consistent with observations of flare-related changes to photospheric 
magnetic fields, which in many cases show fields becoming “more horizontal” ( e.g ., Wang and Liu 2010). We revisit 
this point below. 
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Given the relative noise levels on the normal and horizontal fields for this magnetogram, it could be argued that the 
Dirichlet and combined fields are too inconsistent with the observed normal field, so discrepancies in their field energies 
are unimportant. The noise-weighted potential field, however, is statistically consistent with the observed normal field, 
but the energy in the Neumann potential field is more than 10 32 erg larger. This discrepancy is a significant amount of 
energy, commensurate with that released in large flares that produce fast CMEs ( e.g., Emslie et al. 2012). Evidently, 
the choice of BCs could have significant consequences for modelers attempting to determine free energies in coronal 
fields. 

It should be noted, however, that our use of DFTs can produce potential energies systematically larger than other 
methods (e.g., a Green’s function approach; Sakurai 1982). This is because our AR is implicitly assumed to be 
surrounded by neighboring copies of itself, meaning its field is partially confined by the magnetic pressure from its 
hypothetical neighbors. This can produce a volume-averaged B 2 higher than that for an isolated AR. But this bias 
should not affect the ordering in energies yielded by each approach. 

Based upon Equation (24), it is plausible that the energy differences between the different potential fields might be 
related to differences in the amount of unsigned vertical flux in each potential model. Accordingly, we summed the 
unsigned vertical magnetic flux in each model over all pixels with \B Z \ > 100 Mx cm 2 , which gave 2.3 x 10 22 Mx for 
the Dirichlet field, 2.8 x 10 22 Mx for the Neumann field, and 2.2 x 10 22 Mx for the combined field. This threshold was 
chosen to avoid including \B Z \ from weak-field regions in the Dirichlet solution, caused by spurious values of • Bh 
from differentiating noisy data there. With this threshold, the Dirichlet summation included about 5.3 x 10 4 pixels, 
the Neumann summation included about 3.8 x 10 4 pixels, and the combined solution included about 3.3 x 10 4 pixels. 
The larger number of pixels in the Dirichlet sum probably reflects greater contributions from noise relative to the sums 
from other fields. With this caveat in mind, we note that the ordering of field energies does not follow the ordering 
of fluxes in each model. The differences between the observed B z and the Dirichlet field’s B z shown in the lower-left 
panel of Figure 1 imply that the Dirichlet solution’s vertical field is systematically weaker in regions where the observed 
vertical field is strong. This, we believe, explains why the Dirichlet solution has less unsigned magnetic flux than the 
Neumann solution, despite the I3 z -halos in the former. 

In the top panel of Figure 7, we show, as functions of time, the widths from Gaussian fits to the cores of histograms 
of B z and Wh- Bh, where the core regions are defined to be within ±20 Mx cm' 2 and ±60 Mx cm -2 pix~ 4 , respectively. 
A linear approximation to the median longitude of the active region’s pixels is plotted on the top axis of the bottom 
plot. (The active region’s latitude remained about S20 throughout.) We assume core widths measure uncertainties in 
B z and • Bh, and fitted widths increase toward the limb, as expected. The bottom panel shows magnetic energies 
(left axis) in the Neumann and noise-weighted hybrid potential fields, and their ratio (right axis). Since noise in the 
Neumann condition increases toward the limb, the noise-weighted potential field incorporates more of the Dirichlet 
information toward the limb; and since the Dirichlet energy is systematically lower than the Neumann energy, the 
ratio of Neumann to noise-weighted energies increases. Over most of the interval plotted, the energy of the Neumann 
potential field is 10% (or more) higher than that of the noise-weighted field, implying estimates of free energy would 
be about 10% (or more) lower if the Neumann field is defined as the minimum energy state. Given the variation in 
weights versus disk position in the noise-weighted approach, we favor making comparisons between potential energies 
at different disk positions determined using fixed weights, as opposed to the position-dependent weights here. 

4.2. Do Other Magnetographs Show Similar Features? 

We have identified two notable patterns in Dirichlet potential fields extrapolated using SDO/HMI vector field mea¬ 
surements: opposite-polarity halos around strong-field regions (middle-left panel of Figure 1); and magnetic energies 
substantially lower than Neumann potential fields. 

Since these might arise from systematic effects in the observation and inversion processes used to estimate the mag¬ 
netic field by the SDO/HMI Team, we also analyze magnetograms of two other active regions, observed by the Synoptic 
Optical Long-term Investigations of the Sun/ Vector- SpectroMagnetograph (SOLIS/VSM; Keller, Harvey, and Giampapa 
2003) and by the Solar Optical Telescope /SpectroPolarimeter (SOT-SP) aboard Hinode (Lites et al. 2013). 

AR 11117 was observed by SOLIS/VSM on 27 October 2010. Tadesse et al. (2012) extrapolated the coronal magnetic 
field in this region to determine the free magnetic energy content before a C-class flare near 17:00 UT. Here, we compute 
the pre-flare potential field energy from the VSM magnetogram recorded around 16:33 UT, which is available from the 
National Solar Observatory’s website. We used a Lambert equal-area projection to remap the data onto a Cartesian 
plane, with 1" pixels. In Figure 8, we show the observed and Dirichlet potential fields. As with the SDO/HMI data 
in Figure 1, opposite-polarity halos in B z are present around strong-field regions in the Dirichlet solution. 
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Figure 7. Top: Widths from Gaussian fits to cores of histograms of B z (solid) and Vh • Bu (dashed; cores are ±20 Mx cm -2 
and ±60 Mx cm -2 pix -1 , resp.) versus time for three days of vector SDO/HMI magnetograms. (Vertical axis uses G instead 
of Mx cm -2 for readability.) Bottom: Magnetic energies (left axis) in Neumann (solid) and noise-weighted hybrid (dashed) 
potential fields, and ratio of energies (±’s, right axis). A linear approximation to the median longitude of the active region’s 
pixels is plotted on the top axis of the bottom plot; the active region’s latitude remained about S20 throughout. We assume core 
widths are measures of uncertainty, and widths for B z (top plot, solid) increase toward the limb. The noise-weighted potential 
field incorporates more of the Dirichlet information toward the limb, so the ratio of energies increases. 


Tadcsse et al. (2012) adopt noise levels of 1 Mx cm' 2 and 50 Mx cm -2 for the line-of-sight and transverse field 
components, respectively. Given the active region’s location near N21W25 at the time of these observations (about 
32° from disk center), transverse field measurements contribute about 25% to the radial field estimate. Hence, crude 
estimates of errors in the radial and horizontal components would be about 12.5 Mx cm -2 and 37.5 Mx cm' 2 , 
respectively, for field components in a plane through disk center. (The horizontal error is probably higher, because 
one component of the horizontal field does not include any line-of-sight component. We have also neglected any 
discussion of ambiguity resolution errors.) Assuming a 3:1 weighting in favor of the radial field measurement, we 
find the Neumann, Dirichlet, and noise-weighted fields have energies of 3.3 x 10 32 erg, 8.6 x 10 31 erg, and 2.4 x 10 32 
erg, respectively. The Dirichlet energy is smaller by a factor greater than 3, so the difference between the Neumann 
and noise-weighted energies is nearly 10 32 erg, about half the free energy estimated (Tadcssc' et al. 2012) from their 
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Figure 8. Left: Grayscale shows vertical field B z in AR 11117 observed by SOLIS/VSM on 27 October 2010 near 16:33 UT. 
Saturation is set to ±125 Mx cm' 2 . Arrows (white where B z < 0 and black where B z > 0) show the observed horizontal field. 
Right: Grayscale shows vertical field B z from a potential extrapolation using the Dirichlet BC determined from the measured 
horizontal magnetic field. Arrows (white where B z < 0 and black where B z > 0) show the potential horizontal field. The x and 
y axes are in interpolated, reprojected HMI pixels, which are square and 362 km on a side. Note the opposite-polarity halos 
around strong-field regions in the Dirichlet field, as in the SDO/HMI data. 


non-linear force-free-field (NLFFF) extrapolation. We remark that all of our potenetial field energies are significantly 
smaller than the potential-field energy estimate of adesse et al. (2012), by roughly a factor of 3. This might be due 
to the much larger field-of-view over which they estimated the magnetic energy. 
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Figure 9. Left: Grayscale shows vertical field B z in AR 10930 observed by Hinode/ SOT-SP on 12 December 2006 near 21:00 
UT. Saturation is set to ±750 G. Arrows (white where B z < 0 and black where B z > 0) show the observed horizontal held. 
Right: Grayscale shows vertical held B z from a potential extrapolation using the Dirichlet BC determined from the measured 
horizontal magnetic held. Arrows (white where B z < 0 and black where B z > 0) show the potential horizontal held. Pixels 
are 0.63" on a side. Note the opposite-polarity halos around strong-held regions in the Dirichlet held, as in the SDO/HMI and 
SOLIS/VSM data. 


Schrijver et al. (2008) analyzed NLFFF extrapolations of NOAA AR 10930 based upon a vector magnetogram 
derived from Hinode/ SOT-SP measurements near 21:00 UT on 12 December 2006. In the left panel of Figure 9, we 
show the observed vector field ( B z in grayscale, horizontal field with vectors). (Pixels are 0.63" on a side.) In this 
figure’s right panel, we show the Dirichlet field (again, B z in grayscale, horizontal field with vectors). As with the 
SDO/HMI and SOLIS/VSM data, opposite-polarity halos are seen in the vertical field around strong-field regions in the 
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Dirichlet solution. Perhaps not coincidentally, the large-scale, strong-field polarities in the observed field (left panel) 
are surrounded by the many small-scale concentrations of opposite-polarity flux. It is possible that the photospheric 
magnetic field dynamically creates an opposite polarity moat via a relaxation process toward a field that more closely 
resembles the Dirichlet potential field. The observed opposite-polarity moats are not, however, as intense as those in 
the Dirichlet solution. 

As with the SDO/HMI and SOLIS/VSM observations, the energy of the Neumann potential field for the 
Hinode/ SOT-SP data, ft 1.6 x 10 33 erg, is also larger than the energy of the Dirichlet field, ft 8.1 x 10 32 erg, a 
difference of a factor near two. With a weighting of 10:1 in favor of the normal-field (Neumann) BC, the weighted 
field’s potential energy is ft 1 x 10 32 erg less than the Neumann field’s. This is roughly equivalent to the ft 10 32 erg 
(Emslie et al. 2012) estimated to have been released in an X-class flare and CME a few hours after this magnetogram 
was recorded. Again, we remark that different choices of BCs can significantly affect estimates of free energies in 
coronal fields. 


4.3. Hints about Photospheric Horizontal Currents 

Beyond the quantitative differences between potential fields inferred using Neumann condition and hybrid BCs, a 
qualitative difference in the hybrid approach is its treatment of the observed B z as neither unrelated to the presence 
of currents nor “gospel” data that should be precisely matched when specifying the potential field. 

In the hybrid approach, differences between the observed and potential B, both account for uncertainties in the 
measurements and indicate that horizontal currents are present at the photosphere. We can define the total field as 
the sum of potential and non-potential parts, 

B tot = B P + jgNP _ (25) 

Then B^ p is related to horizontal current, 

4t xJ x /c=d y B z - d z B v = d v B™ p - d z B" p (26) 

4tt J y /c= d z B x - d x B z = d z B™ p - d x B* p . (27) 

Magnetic field inversions that estimate the field at a single height cannot provide information about variation in 
magnetic field components along the LOS (or z) directions, so the structure of horizontal currents cannot be determined 
from these observations alone. 

Nonetheless, patterns in the spatial structure of B^ p might reveal clues about horizontal photospheric currents. 
For instance, we have noted the opposite-polarity halos surrounding strong-field regions in the hybrid (or Dirichlet) 
potential fields. When differencing the observed and hybrid (or Dirichlet) potential normal magnetic fields in a sunspot 

as in, for instance, the bottom-left panels of Figures 1 and 2 — the halos around the positive-polarity spot at the 
upper right of the magnetogram produce an interesting structure for B^ p : a core field with opposite polarity to that 
of the spot, surrounded by a moat of the same polarity as the spot. This structure would be produced by a horizontal 
electric current flowing clockwise (CW) around the spot’s periphery. This current would flow in the opposite sense 

of the counterclockwise (CCW) sheath current expected to flow around an isolated, positive flux tube surrounded by 

field-free plasma ( e.g ., Spruit 1981). Bommier et al. (2011) have also reported similar current structures in another 
active region, deduced using a different approach: multiple-height vector magnetic field measurements. 

It should be noted that our interpretation of these non-potential field patterns as manifestations of horizontal 
currents implies that they are distinct from non-potential field patterns associated with vertical currents observed to 
cross the photosphere (e.g., sheared, strong-field polarity inversion lines; Georgoulis, Titov, and Mikic 2012). Hence, 
the relationship between the horizontal currents that we describe and vertical currents, whether outgoing/ “direct” or 
inflowing/“return” currents (e.g., Torok et al. 2014), is unclear. 

It is also possible that the halos arise from the unphysical nature of the potential field assumption: the actual 
active-region field is confined by surrounding, quiet-sun regions, but potential fields are by definition force-free, so 
are space-filling. In fact, Neumann potential fields typically exhibit this effect: Neumann solutions tend to have 
strong, spatially coherent horizontal magnetic fields outside of intense concentrations of vertical field. This contrasts 
to the observed horizontal Helds, which tend to be much more strongly confined to areas with strong vertical flux: the 
observed horizontal fields in Figure 10, where \Bh\ from the SDO/HMI dataset is shown in grayscale, closely match 
areas where strong vertical fields are present in the upper-left panel of Figure 1. A moat of extended horizontal fields 
from the Neumann solution can also be seen in Figure 10, where the 200 Mx cm -2 contour line of the Neumann field’s 
\Bh\ (the outermost red contour) lies well outside the same-strength contours of both \Bh\ observed by SDO/HMI 
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(aqua contour) and the Dirichlet field | Bh\ (blue contour). These “excess” horizontal fields in the Neumann solution 
indicate that, relative to the potential-field model, the actual magnetic field is “too vertical.” The Dirichlet field’s 
dearth of vertical flux (compared to the observed field) accords with this characterization: relative to the vertical field 
in the Dirichlet model, which is derived from the observed horizontal field, the actual field’s area-integrated \B Z \ is 
higher, i.e., the real field tends to be more vertical. 



Figure 10. The grayscale shows the smoothed horizontal field strength \Bh\ from the SDO/HMI observation in Figure . The 
x and y axes are in interpolated, reprojected HMI pixels, which are square and 362 km on a side. The aqua contour line shows 
the 200 Mx cm -2 level curve of the observed \Bh\- The red and blue contour lines show 200 and 400 Mx cm' 2 level curves of 
\Bh\ from the Neumann and Dirichlet potential fields, respectively. Saturation of the grayscale is set to 400 Mx cm -2 . Note 
the presence of significant horizontal fields in the Neumann solution outside of the observed horizontal field (i.e., outside the 
aqua line). 


We suspect that opposite-polarity the halos present in the Dirichlet field’s (modeled) B z likely arise from the rapid 
fall-off observed in horizontal field strengths, which tend to produce halos in V/, • Bh like those discernible in the 
middle-right panel of Figure 1. Simply put: halos in the observed • Bh produce halos in the Dirichlet model’s B z . 

Despite the presence of opposite-polarity halos elsewhere in Figures 1-2 and in Figures 8-9, similar structures in 
B^ p are not ubiquitous in other areas of the magnetograms we have analyzed here. Studying Dirichlet / hybrid poten¬ 
tial magnetic fields in additional active regions, as well as additional multi-height vector magnetic field measurements, 
would be useful to better understand the structure of electric currents within active region fields. 

5. SUMMARY AND DISCUSSION 

We have presented methods to compute potential (current-free) magnetic fields in the photosphere and corona that 
incorporate measurements of the horizontal photospheric magnetic field. The horizontal photospheric field determines 
a 2D Dirichlet boundary condition at the photosphere to solve Laplace’s equation in 3D. In addition to matching just 
the horizontal field, we also described methods that attempt to match both vertical and horizontal components of the 
field, in a statistical sense. 

Historically, potential magnetic fields have been computed using estimates of just the vertical photospheric magnetic 
field, which determines a Neumann boundary condition there for Laplace’s equation in three dimensions. This was 
perhaps because the line-of-sight magnetic field can be measured more easily and accurately than the component of 
the magnetic field transverse to the LOS, and near disk center the LOS field approximates the radial field. In fact, 
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the difficulty of measuring the full magnetic vector at the photosphere meant that such observations were relatively 
rare prior to the advent of SOLIS/VSM, Hinode/ SOT-SP, and SDO/HMI. Given the expected presence of horizontal 
electric currents in active regions at the photosphere (indeed, throughout the photosphere-to-corona volume), however, 
there is no physical basis for assuming that a potential field’s vertical component should precisely match that of the 
observed field. Hence, the customary use of Neumann boundary conditions for potential field extrapolations likely 
arose from observational capabilities, not physical considerations. 

We showed that a hybrid field, formed by a linear combination of separate solutions to Laplace’s equation using 
the Neumann and Dirichlet boundary conditions, matches the observed vector magnetic field more closely than either 
the Neumann or Dirichlet solutions alone. If the aim of determining a potential field is to find the current-free coronal 
magnetic field most consistent with the observed photospheric field, then evidently information about the horizontal 
field should also be included when deriving the potential field. Differences between the observed field and the hybrid 
field that most closely matches it can arise from noise in the magnetic field measurements, but are also expected in 
the presence of electric currents. 

We found two systematic features of Dirichlet and hybrid fields. First, we noted that opposite-polarity halos of 
vertical magnetic field surrounding strong-field concentrations of vertical flux are typical in the magnetograms we 
studied. These are not seen in the observed vertical field. Second, we found that the magnetic energy of the Dirichlet 
and hybrid potential fields tends to be significantly lower than that of the Neumann potential fields. The differences 
in magnetic energies of these potential fields - 10 32 erg or more — are commensurate with estimates of the energy 

released in solar eruptive events ( e.g ., Ernslie et al. 2012). (As noted above, the assumption of periodicity implicit in 
our use of DFTs to compute potential fields might yield energies that are artificially large.) We found that Dirichlet or 
hybrid potential fields derived from SDO/HMI, SOLIS/VSM, and Hinode/ SOT-SP vector magnetograms all exhibited 
both properties. 

Our results have two notable implications for estimates of free magnetic energy in non-potential coronal magnetic 
field models. First, from a physical point of view, the free magnetic energy necessarily depends upon the assumed 
minimum-energy field, and hybrid potential fields both match observations more closely than Neumann potential 
fields and have lower energy. This suggests that estimates of coronal free energy made using the Neumann field as 
minimum-energy field might understate the true free energy significantly. 

Second, from a practical standpoint, hybrid potential fields probably provide better initial conditions than Neumann 
fields for some methods to calculate model coronal fields. For instance, the initial coronal field in the optimization 
method {e.g., Wheatland, Sturrock, and Roumeliotis 2000; Wiegelmann 2008) of deriving NLFFFs is typically extrap¬ 
olated from a photospheric Neumann condition, but the relaxation is initiated by replacing the potential photospheric 
boundary field with the observed field (or a “preprocessed” version of the observed field). A large discrepancy between 
the observed boundary field and the model field’s first interior layer is then present. The optimization algorithm 
then attempts to remove any Lorentz forces or magnetic divergences due to this discrepancy by iteratively altering 
the model field above the photosphere. Because hybrid-potential photospheric fields more closely match the observed 
photospheric field, the optimization algorithm should converge more quickly starting from a hybrid potential solution, 
since discrepancies between the initial and observed V/, • B h should be smaller. This can be understood by consid¬ 
ering a Helmholtz decomposition of Bh, which expresses Bh in terms of a scalar potential and a stream function: 
the scalar potential is derived from ■ Bh, and can be used to determine the irrotational part of Bh', the stream 
function is derived from the solenoidal component of Bh, from z ■ (V/j x Bh), and can be used to determine the 
remaining part of Bh (which is directly related to the vertical photospheric electric current). (This decomposition 
underlies a poloidal-toroidal decomposition [PTD] of the full, three-component vector, B; Fisher et al. (2010) and 
(Kazachenko, Fisher, and Welsch 2014) discuss PTD at length, and include detailed descriptions of methods to de¬ 
termine the scalar potential and stream function.) Essentially, NLFFF methods should be focused on matching the 
observed vertical currents, related to z ■ (V/, x B h ), not discrepancies between a Neumann-derived V/, • B h versus the 
observed V/,, • Bh- In fact, it is likely that discrepancies in V/j, • Bh introduce significant problems in NLFFF extrapo¬ 
lations: Valori et al. (2013) report that errors in satisfying V • B = 0 in extrapolated fields can compromise magnetic 
energy estimates, and starting from a potential extrapolation that is more consistent with the observed \7 h • Bh 
should minimize introduction of spurious divergences. For similar reasons, we expect that using hybrid-field initial 
states in magnetofrictional methods {e.g., Valori et al. 2010) could also result in closer agreement between observed 
and extrapolated fields. 

As another consequence of using a hybrid potential field, some treatments of magnetic helicity, such as those 
developed by Berger and Field (1984) and and Antonsen (1985), are not applicable with hybrid fields. These 
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treatments sought to define a gauge-invariant hclicity for a given magnetic field, B, by defining that field’s hclicity 
“relative” to the “reference” potential magnetic field satisfying the corresponding Neumann boundary condition. For 
these definitions of hclicity, precise agreement between the actual and potential fields’ normal components is required 
for gauge invariance. We note, however, that other formalisms that yield well-defined values for magnetic hclicity have 
been developed ( e.g ., Hornig 2006; Low 2015) that do not require any reference field. 

For simplicity, we used DFTs to compute potential fields discussed above. In the Appendices, we present methods for 
the more general Cartesian cases of non-periodic Neumann and Dirichlet potential fields, using finite-element sources. 
For the global spherical case, Nemenman and Silbergleit (1999) present Green’s functions for the 3D Neumann and 
Dirichlet boundary conditions (without a source-surface boundary condition), which can be used to match observations 
with arbitrary spatial resolution. To derive the Dirichlet condition from observations, the 2D Poisson problem in 
spherical coordinates can be solved using FISHPACK, a suite of Fortran codes developed at NCAR by Schwarztrauber 
(1975). An IDL wrapper for FISHPACK has been developed as part of the CGEM project (http://cgem.stanford.edu/), 
and should be publicly released soon. We note that, in principle, spherical harmonic transforms could be used for 
either of these 3D or 2D problems, but the numerical expense of computing very-high-order expansions to match 
high-resolution observations can be prohibitive. 
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APPENDIX 

A. DIRICHLET SOLUTION FOR A TILE 


For some applications, the periodicity inherent in using DFTs to calculate 3D potential fields from the 2D Dirichlet 
BC is problematic. Here, we outline an approach to determine a non-periodic, 3D potential function in Cartesian 
geometry. In a domain that is finite in x and y , one can calculate the 2D Dirichlet BC, x{x, y , 0), for the 3D potential 
function x(x, Vi z ) directly from the boundary values of V/j • Bh(xi,yi) | z =o, without solving the 2D Poisson’s equation. 
The approach is similar to using a Green’s function: superposition is employed with the fundamental solution to the 
2D Laplace equation, 

XD,pt (x,y) = ln(v /(x - x') 2 + (y- y') 2 ) , (Al) 


which gives the potential at (x,y) due to a point divergence of strength D pt at (. x',y'). The divergence of the 2D 
vector field computed from the horizontal gradient of this potential vanishes for all points other than [x',y'). With 
the horizontal divergence expressed as a 2D function D(x,y) = V/ t • Bh(xi,yi)\ z —o, the potential at z = 0, xv(x,y,0) 
is found by summing the contributions from all sources over the 2D domain. (Unlike the situation with a Green’s 
function, an integral over the fundamental solution does not converge on an infinite domain.) 

Once the potential function is specified on z = 0, the potential at an arbitrary point in the domain (x, y, z ) can be 
computed via 


X{x, y, z ) 


tJ 


tL4' 


x( x ',y', o) 


[{x — x') 2 + (y — y') 2 + 2 2 ] 3 / 2 


(A2) 


It should be noted that the potential at given height H above a point (xg,yo) on the surface will be significantly 
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influenced by values of the potential on the surface at least a horizontal distance H away from (xo,yo). Therefore, to 
extrapolate the field above a rectangular region of the photosphere N x x N y pixels in extent to a height commensurate 
with the area, H ss N x x N y , the potential xd(^, y , 0) must be known over a substantially larger area of the z = 0 
plane — of order 3 N x x 3 N y or more. In contrast, when extrapolating from a Neumann BC in a region, the normal 
field exterior to that region is often assumed to vanish, meaning the Neumann extrapolation does not require input 
from so large an area. Hence, the Dirichlet field’s dependence upon the potential at the surface, which in general 
extends far outside the Neumann field’s sources, is a disadvantage in practical terms. One possible work-around is to: 
determine xd(£ ; y , 0) over an area only slightly larger than the N x x N y region; extrapolate upward just enough zones 
to determine d z \{x , y , 0) from a one-sided, finite-difference stencil for the vertical derivative; and then use this normal 
derivative (not the observed normal field) as a Neumann boundary condition to extrapolate the potential field. 

One practical problem in computing xd(£, 2/,0), however, is that the fundamental solution is singular at a given 
source’s location. For numerical determination of Dirichlet potential fields (for instance, to integrate field lines), this 
is a major shortcoming. To overcome this deficiency, one can treat the horizontal divergence computed for the pixel 
Dij at (xi,yj) as that due to an extended source with spatially uniform areal density D = Dij/(As) 2 over a square 
“tile” of finite extent As in x and y, corresponding to the pixel area. The resulting expression for the potential due to 
this finite source is not singular at the pixel center, or elsewhere (for a finite domain in x and y). 

We now determine the finite-element potential, assuming the pixel is centered at the origin. It suffices to compute 
the potential at two points: (i) at the center of the tile; and (ii) at an arbitrary point outside the tile, which we will 
use to compute the potential at neighboring pixels. For points far from the tile, relative to As, we can approximate 
the potential as that from a point source. 


A.l. Solution at Tile Center 

We first want to compute the potential corresponding to the point at the center of the tile, 

r As/2 r As/2 


/ LAS/ A /•LAS/ Z _ 

dx' / dy' —— ln(\/ x a + y' 2 ) , 

-As/2 J — As/2 


(A3) 

s/2 J — As/2 

where the argument of xd refers to distance with respect to the center of the source pixel, the origin for x' and y'. 
Then 

jj rAs/2 fAs/2 


_ ta pLAS/Z pLAS/Z 
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—^ {2 ln(As) - ln(2) - 3 + tt/2} = {2 ln(As) - 2.122} 

47r 47r 


(A12) 


A.2. Solution Outside of Tile 

For the potential exterior to the tile, but near it, we now want to compute 

Xd (x,y)= [ dx' j dy'^- ln(y/( x - x') 2 + (y - y') 2 ) , (A13) 

J —As/2 J~As/2 Z7T 

where, as above, the arguments of xd refers to distance with respect to the center of the source pixel, which we take 
as the origin. 

We first change variables, x" = x — x' and y" = y — y' , and define 


x± = x ± As/2 
y± = y± As/2 . 


(A 14) 
(A15) 


Then the upper and lower limits of integration ±A/2 for x' become x=p for x", and similarly for y' and y". We then 
have 


Xd (x,y) = —^ J dx" J dy" ln(x" 2 + y" 2 ) 
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Xd (x,y) = 
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as 


This expression can be recast in terms of the measured horizontal divergence Dij at (xj, yj) and pixel indices (i, j) 
Xd (x,y) = (2 ln(As) - 2 
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Figure Al. Grayscale map of B z in AR 11158 from the Dirichlet BC computed using the finite-element (“tile”) approach 
derived in Appendix . The x and y axes are in interpolated, reprojected HMI pixels, which are square and 362 km on a side. 
Saturation is set to ±1000 Mx cm - . The opposite-polarity halos present around strong-field regions in the Dirichlet magnetic 
field computed with DFTs (see middle-left panel of figure 1) are also visible here. 


We give an example of B z computed with the finite-element approach in Figure Al, which shows the Dirichlet field’s 
B z in AR 11158. The method is numerically expensive (compute time scales as iV 4 for each 2D slice in z) : so we 
only computed the held on a 426 x 276 cropped region of the full 621 pixel x 610 array. We exclude a 10-pixel-wide 
strip around the cropped array’s perimeter from this image, since the solution in much of this strip was strongly 
influenced by edge effects (owing to the cropping). This solution also exhibits the opposite-polarity halos present 
around strong-held regions that were visible in the Dirichlet held computed with DFTs (see middle-left panel of Figure 
!)• 


B. NEUMANN SOLUTION FOR A TILE 


Using hybrid Neumann-Dirichlet potential helds requires computing the Neumann and Dirichlet potentials with 
consistent approaches. Hence, in addition to the method presented above for computing the Dirichlet potential from 
a hnite source — a uniform-density tile — we also describe a method to compute the potential arising from a tile of 
vertical flux, the analogous source for the Neumann case. As above, we use superposition to compute the potential 
from tile, building upon the Green’s function for the 3D Laplace equation, 


XN.pt (x, y,z) 


_L_<hpt_ 

2tt yj(x — x') 2 + {y — y') 2 + z 2 


(B23) 


which gives the potential at (a ’,y,z) due to a point source of hux 4> p t at 0). As with the Dirichlet case, this 

potential is singular at a point source’s location, which can be problematic for numerical determination of Neumann 
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potential fields. We again treat the flux in a pixel at (xi,yj) as that due to an extended source with spatially 
uniform areal density B z = <Fy/(As) 2 over a square tile of finite extent As in x and y, corresponding to the pixel 
area. As with the Dirichlet case, the resulting expression for the potential due to this finite source is not singular at 
the pixel center. 

We now determine the finite-element potential, assuming the pixel is centered at the origin. It suffices to compute 
the potential at two points: (i) at the center of the tile; and (ii) at an arbitrary point outside the tile, which we will 
use to compute the potential at neighboring pixels. For points far from the tile, relative to As, we can approximate 
the potential as that from a point source. The total potential arising from an arbitrary distribution of fluxes over a 
set of pixels can then be computed by summing the potentials from all tiles. 


B.l. Solution at Tile Center 

We first want to compute the potential corresponding to points at the center of the tile, 


o pAs/2 rLAs/'Z 
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This result implies that, to compute the potential at the point corresponding to the center of a pixel, assuming a 
uniform distribution of flux over that pixel, the denominator in Equation (1123) — which is zero at that point — 
should be replaced with As/(41n[-\/2 + 1]). 


B.2. Solution Outside of Tile 

For the potential exterior to the tile, but near it, we now want to compute 


B~ f As / 2 f‘As/2 1 

Xn (x.y.z) = —A / dx' / dy '— . = , 

2tt J-as /2 J-As /2 y/{x - x') 2 + {y - y') 2 + z 2 


(B31) 


where, as above, the arguments of xn refers to distance with respect to the center of the source pixel, which we take 
as the origin. We change variables via Equations (Al l) and (A15), so the upper and lower limits of integration ±A/2 
for x' become for x", and similarly for y' and y" 1 so 


XN(x,y) 
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Defining a = ( x" 2 + z 2 ), 
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(B33) 
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(B35) 


We had no success at further reducing this expression analytically, so instead evaluated it numerically for points in 
the neighborhood of the tile. For five classes of nearby, off-tile points, discrepancies between the point-source potential, 
Equation (B23), and our numerical results were above 1.5%. For these sets of points, we used our results to derive 
coefficients to “correct” the distance factor in the denominator of Equation (B23). For two points directly above the 
tile’s center, at z = 1 and z = 2, the denominator should be multiplied by 0.9286 and 0.9800, respectively. For the 
four nearest neighbors in x and y, the denominator should be multiplied by 1.0381 for z = 0 and 0.9878 for z = 1. 
Finally, for the four corner neighbors in x and y at z = 0, the denominator should be multiplied by 1.0249. 

Note that using these “corrections” for 0 < z < 2 will produce a Neumann potential that disagrees in these zones 
with the Dirichlet potential, computed via Equation (A2), extrapolated from the Neumann potential at z = 0. 


1 By differentiating Equation (B 3) prior to integrating the source distribution for a tile, Barnes et al. (2009) were able to derive an 
analytic expressions for the field’s components; see http://www.cora.nwra.com/~graham/Potential/ 
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